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THE DETAILED FUNDAMENTAL PLANE OF TWO HIGH REDSHIFT CLUSTERS: MS 2053-04 AT Z=0.58 
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ABSTRACT 

We study the fundamental plane relation in high redshift clusters using a sample of 26 galaxies in 
MS 2053-04 (z = 0.583) and 22 galaxies in MS 1054-03 (z = 0.83). The zeropoint and scatter are 
compared to results for lower redshift clusters in order to trace evolutionary effects. Furthermore, our 
large sample enables us to investigate correlations between residuals from the fundamental plane and 
other characteristics of the galaxies, such as color, HP linestrength, spatial distribution, and mass. 
The observed scatter of the early-type galaxies with a > 100 km s~ 4 around the fundamental plane 
is 0.134 and 0.106 in logr e for MS 2053-04 and MS 1054-03 respectively. The residuals from the 
fundamental plane of MS 2053-04 are correlated with residuals from the HP — a relation, suggesting 
that stellar populations are playing a role in shaping the fundamental plane. The measured evolution 
in log M/L is influenced by selection effects, as galaxies with lower M/L in the Johnson B-band enter 
a magnitude-limited sample more easily. When we select high mass early-type galaxies to avoid this 
bias, we find logAf/L^ ~ — 0.47z and a formation redshift Zf orm ~ 2.95, similar to earlier results. 

Subject headings: galaxies: clusters: individual (MS 2053-04 , MS 1054-03 ) 
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1. INTRODUCTION 

In the local universe early-type galaxies follow a tight 
scaling relation, known as the Fundamental Plane (FP). 
The relation between effective radius, central velocity 
dispersion and surface brightness r e ~ o- a I@, which is 
a plane in (log r e , log er, log/ e ) space, was discovered by 
Djorgovski & Davis (1987) and Dressier et al. (1987). In 
combination with the virial theorem 



M/L 



o~ 2 r 



(1) 



the small scatter around the FP implies that, under the 
assumption of homology, the M/L ratios of early- type 
galaxies are well behaved and scale as 



M/L 



(2) 



As the M/L ratio increases with ageing of a stellar pop- 
ulation, the FP is a useful tool in research on galaxy for- 
mation and evolution. Based on a sample of 226 E and 
SO galaxies in 10 clusters of galaxies, J0rgensen, Franx, & 
Kjaergaard (1996, hereafter JFK96) concluded the local 
plane in the Johnson B band has the form: 



logr e = 1.201ogcr- 0.83 log I e +7 



which implies that 



M/L ~ M 



0.28 r -0.07 



(3) 



(4) 



Later studies on intermediate and high redshift clusters 
of galaxies used the zeropoint shift of the plane to es- 
timate average formation redshifts of the stars in early- 
type galaxies (e.g., Bender et al. 1998; van Dokkum et al. 
1998, hereafter vD98; J0rgensen et al. 1999; Kelson et al. 
2000c, hereafter K2000; Pahre et al. 2001; van Dokkum 
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& Stanford 2003). The scatter around the plane provides 
constraints on the spread in galaxy ages. 

The slope of the FP (and other scaling relations) con- 
strains systematic age trends with mass and other pa- 
rameters. Any evolution of the slope of the FP with red- 
shift implies that the ages of the stellar populations are 
correlated with galaxy mass. For the cluster CL1358+62 
at z = 0.33 K2000 finds that the slope of the FP has 
not evolved significantly over the past ~ 4 Gyr. The 
sample of 5 bright MS 2053-04 galaxies used by Kel- 
son et al. (1997, hereafter K97) to study the FP at z 
= 0.583 seemed to agree with the values of coefficients 
a and P as given by JFK96. However, the sample was 
too small to perform a proper fit. The same conclu- 
sions were drawn for MS 1054-03 at z = 0.83 based on 
6 early- type galaxies (vD98). In this paper we investi- 
gate the FP of MS 2053-04 and MS 1054-03 early-type 
galaxies using larger samples spread over a larger range 
of distances from the brightest cluster galaxies (BCG) . In 
Sect.|5]we discuss the spectroscopy, sample selection and 
velocity dispersions. Imaging and measurement of the 
structural parameters is described in Sect. Zeropoint 
of the FP with JFK96 coefficients and scatter around the 
plane are discussed in Sect. QJ In Sect. [S] we study cor- 
relations between the residuals from the FP and various 
other properties of the galaxies. Finally the conclusions 
are summarized in Sect. VEGA magnitudes are used 
throughout this paper. We use Ho = 70 km s~ x Mpc -1 , 
Q\ = 0.7, and Qm = 0.3, but note that our results are 
independent of the value of the Hubble constant. 



2. SPECTROSCOPY 



2.1. 



Sample selection and observations 

All spectra used in the FP analysis were obtained with 
the LRIS spectrograph (Oke et al. 1995) on the 10 m 
W.M. Keck Telescope. The data were obtained in a 
series of observing runs from 1996 to 2002. The ma- 
jority of galaxies in our final sample were selected on 
the basis of their spectroscopic redshift and their I- or 
F814W-band magnitude. The samples were limited at 
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/ w 22 for all runs; galaxies with 7 < 21 were given 
highest priority in the mask designs. The redshift infor- 
mation came from a large /-selected spectroscopic sur- 
vey of both clusters described in detail in Tran et al. 
(1999), van Dokkum et al. (2000), and Tran (2002). For 
the initial MS 1054-03 observing runs only limited red- 
shift information was available, and we applied color cri- 
teria to select likely cluster members. Galaxies having 
A(R -I)- 0.25A(S -R) < -0.4, with A(R - I) and 
A(B — R) colors relative to the central galaxy, were ex- 
cluded. The color ranges were chosen such that blue 
cluster members were unlikely to be excluded. The final 
FP sample of MS 1054-03 contains 19 galaxies that were 
selected with these mild color constraints. No morpho- 
logical information was used in the selection process. 

For most observations we used the 600 lines mm -1 
grating blazed at 7500 A; some of our earlier data were 
taken with the 831 lines mm -1 grating blazed at 8200 A 
(see vD98). The wavelength coverage was typically 
~ 3500 to ~ 5400A in the restframe. Exposure times 
ranged from 7500 s to 33400 s and from 10500 s to 22800 s 
for the MS 2053-04 and MS 1054-03 galaxies respec- 
tively. The instrumental resolution was typically <Ji na t ~ 
40 — 80 km s _1 and signal-to-noise ratios ranged from 20 
to 100A _1 in the observed frame (in the continuum). 

A total of 43 galaxies (26 early-type) were observed 
in MS 2053-04. The MS 1054-03 sample contained 30 
galaxies (14 early-type). The morphological classification 
is described in Sect. 13.31 Early-type galaxies include E, 
E/S0 and SO morphologies. 

2.2. Basic reduction 

The spectra were reduced using our own software and 
standard IRAF software routines (see, e.g. Kelson et 
al. 2000b). The wavelength calibration was performed 
using the night sky emission lines. The typical rms scat- 
ter about the fitted dispersion solutions is about 1/15 
of a pixel. Since the dispersion is about 1.28 A pixel -1 
for the data taken with the 600 mm -1 grating (and 
0.92 A pixel' 1 for the 831 mm -1 data), the rms scatter 
is equivalent to velocity errors smaller than 5 km s _1 . 

The flat-fielding accuracy is generally better than a 
percent, on small scales. The data have not been ac- 
curately flux-calibrated so on large scales, the notion of 
flat-fielding accuracy is not meaningful. There tend to 
be moderate-scale (k = 100 A -1 ) residuals in the flat- 
fielding that are the result of a mismatch between the 
fringing in the flat-fields and the fringing in the data. 
Such inaccuracies in de-fringing the data have no effect 
on the velocity dispersions because the spectra are ef- 
fectively filtered on those scales (and larger) in the pro- 
cess of matching the continua of the template and galaxy 
spectra. 

The subtraction of the sky was performed using stan- 
dard, published, and well-tested methods: for each 
galaxy the two-dimensional spectra were first rectified, 
and then 1st- or 2nd-order polynomials were fit to the 
pixels on both sides of the galaxy, typically excluding the 
few arcsec where the galaxy is bright. An iterative, clip- 
ping routine was used to reject any remaining bad pixels 
from the fit, similar to what IRAF allows the user to do. 
These methods have been discussed elsewhere at great 
length and we choose not to bore the reader with very 



familiar territory. While more complicated and precise 
means of sky subtraction are now available (see Kelson 
2003), these data were analyzed before such methods be- 
came available, and the accuracy of the sky subtraction 
performed in the " traditional" way is satisfactory for our 
purposes. 

The slit widths varied slightly with each run, rang- 
ing from 0'.'90 to l'/05. The extraction apertures for 
MS 2053-04 and MS 1054-03 were 7 CCD rows (or l'/5). 

2.3. Velocity dispersions 

Velocity dispersions were measured with a direct fit- 
ting method (Kelson et al. 2000b). Details of the fit- 
ting procedure are given in Kelson et al. (2003). Direct 
fitting methods provide several advantages over Fourier- 
based techniques. Most importantly, pixels are no longer 
weighted uniformly in the computation of the fitting met- 
ric (in this case \ 2 )- As in Kelson et al. 2000b, we weight 
each pixel by the inverse of the expected noise (pho- 
ton and electronic read noise). Furthermore, the spectra 
that exhibit strong Balmer absorption had those features 
given zero weight in the fitting, in order to minimize the 
contribution of those features to template mismatch er- 
ror, and also to ensure that our estimates of a reflected 
the old, underlying stellar populations. We typically fit 
the spectra over a ~ 1000 A wavelength range (in the 
restframe). We used a range of template stars, from G5 
through K3 and adopted the template star that gave the 
lowest mean \ 2 ■ In both clusters, HD102494, a G9IV 
star, was the "best" template. 

Both observations and simulations by J0rgcnscn et al. 
(1995b) showed that, at low S/N, measured velocity dis- 
persions were systematically too large. At the same S/N, 
this systematic effect was largest for the galaxies with ve- 
locity dispersions below 100 km s _1 . Therefore, we omit 
all galaxies with a < 100 km s _1 from our samples. Fur- 
thermore, we limit our samples to sources with an error 
in a smaller than 15 %. Errors in the dispersion and ve- 
locity are initially determined from the local topology of 
the x 2 (c, V) surface. 

Seven MS 2053-04 galaxies with a > 100 km s _1 were 
observed both during the 1997 and 2001 observing run. 
A direct comparison between the derived velocity dis- 
persions (prior to the aperture correction) is presented 
in Fig. ^ Generally, the agreement is good, with one 
outlier. All spectra have similar S/N. The rms value 
of "•97-0-01 j s i 23, slightly higher than the expected 

V do 97+ <J<T 01 

value of 1. The mean deviation between the two runs is 
— 2 ± 5%, consistent with zero. We conclude there is no 
evidence for systematic effects and conservatively multi- 
ply all errors by the factor 1.23. 

The final sample consists of 26 galaxies (19 early- 
type) in MS 2053-04 and 22 galaxies (12 early-type) in 
MS 1054-03. 

We applied an aperture correction to a nominal aper- 
ture of D = 3'.'4 at the distance of Coma (see JFK96). 
The final velocity dispersions have therefore been multi- 
plied by a factor of 1.057 for MS 2053-04 and 1.062 for 
MS 1054-03. This correction allows for a fair compari- 
son between clusters at a range of redshifts. The data 
are tabulated in Table 

3. IMAGING 
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Fig. 1. — A direct comparison of velocity dispersions for 7 
galaxies in MS 2053-04. Measured a values, prior to aperture 
correction, for the 2001 run are plotted against a values from 
the 1997 spectra. Formal errors derived from X (<x, V) are 
drawn. Based on this overlap sample, actual error bars are 
estimated to be 23% larger. 

For both MS 2053-04 and MS 1054-03 large HST 
WFPC2 mosaics were obtained, each consisting of 6 
pointings. Both clusters were observed in the F606W and 
F814W filters. The layout of the MS 2053-04 mosaic is 
described in Hoekstra et al. (2002). Exposure times were 
3300 s in F606W and 3200 s in F814W per pointing. The 
MS 1054-03 mosaic is described in van Dokkum et al. 
(2000); exposure times were 6500s for each pointing and 
in each filter. Interlacing of the images improved the 
sampling by a factor \/2 for MS 1054-03. 

3.1. Structural parameters 

In this section, we describe the measurement of effec- 
tive radii r e and surface brightnesses at those radii I e . 
For the Johnson B passband, I e in Lq pc~ 2 is related to 
/i e in mag arcsec" 2 as 

log/ e = -0.4(/i e -27.0). (5) 

We created postage stamps sized 12'.'8 x 12'.'8 for the 
26 MS 2053-04 and 22 MS 1054-03 galaxies and fit 2D 
r 1 /™(n = 1,2,3,4) law profiles, convolved with Point 
Spread Functions (PSF), to the galaxy images. As PSFs 
depend on the positions of objects on the CCDs, we used 
Tiny Tim v6.0 to create an appropriate PSF for each 
galaxy. Other parameters determining the shape of the 
PSF are template spectrum (M type star), PSF size (3"), 
sampling and filter (F814W). The code allows simulta- 
neous fitting of the object of interest and any neighbour- 
ing objects. The fits were restricted to radii of 3" to 
5" around the objects, depending on their size. Image 
defects were masked in the fit, as well as neighbouring 
galaxies not well fitted by r 1 / 4 laws. All other pixels got 
uniform weight. 

We performed the r 1 ' n fits for Sersic numbers n = 
1 (exponential), 2, 3 and 4 (de Vaucouleurs law). In 
this paper we always use r e and I e based on r 1 / 4 fits to 
all galaxies on the postage stamps, even if other Sersic 
numbers result in a better \ 2 of the fit. Fitting a r 1 / 4 
profile resulted in a \ 2 < 1-5 for 69% of all galaxies; 86% 
have x 2 < 2. The galaxies in our final samples together 
with the residuals after profile fitting arc presented in 
Fig. H (MS 2053-04) and Fig. [ft MS 1054-03). 



3.2. Error in the structural parameters 

For MS 1054-03 each of the 6 pointings was observed 
twice, with a shift of 0.5 pixels, providing a direct way 
to measure the error in the structural parameters. From 
fits to the two independent observations we infer that 
errors are small (< 4% in r e ) for galaxies larger than 
0'.'83. The error in r e rises to at most 18.5% for smaller 
sources. Though this might seem problematic, we note 
that this large error does not enter the FP analysis as 
the combination r e I~@ enters the FP. Using j3 = —0.83 
from JFK96, the error in the FP parameter r e I®- 83 is 
limited to ~ 2.1% rms, ignoring one outlier with 30% 
offset. The fit of a de Vaucouleurs profile to the out- 
lier ID5347 with merger/peculiar morphology is clearly 
unstable. The uncertainty in the FP parameter is com- 
parable to the ~ 2.5% rms error estimate from Kelson et 
al. (2000a). The small error in the FP parameter is an 
artifact of the slope of the de Vaucouleurs growth curve. 
Hereafter, we use the average of the independent r e mea- 
surements (which reduces the error by l/y/2) and the I e 
corresponding to this average value, using the empirical 
result that r e I® 66 is the most stable combination of the 
structural parameters. In general we can conclude from 
the error analysis that random errors are small. There- 
fore the scatter found around the FP will not be due to 
errors in the photometry. 

3.3. Visual and quantitative classifications 

Galaxies were visually classified by P.G. van Dokkum, 
M. Franx,& D. Fabricant using the procedure as de- 
scribed in Fabricant, Franx,& van Dokkum (2000). In 
this paper we consider cluster members with early-type 
morphology (E, E/S0, SO). Late- type morphologies are 
also plotted, but are not included in fitting procedures 
unless mentionned otherwise. A quantitative alternative 
to the classification by eye is based on the Sersic number 
that results in the smallest \ 2 of the profile fit. 

Fig. 0] shows histograms of the visual classifications for 
the galaxies with best fitting Sersic number 1, 2, 3 and 4. 
As could be expected, early-type morphologies generally 
correspond to Sersic numbers 3 and 4. We note that the 
difference in \ 2 between n = 3 and 4 is often too small to 
consider them as a seperate class of objects. We conclude 
that the visual and quantitative classifications divide the 
sample in roughly the same bulge and disk dominated 
classes. 

3.4. Transformation to restframe magnitude 

In order to compare the FP of MS 2053-04 and 
MS 1054-03 with the FP of Coma at z = 0.023 we 
have to transform the effective radius to units of kpc. 
Furthermore, for a meaningful comparison it is neces- 
sary to study all data in a common photometric band 
in the restframe of the galaxies. The observed F606W 
and F814W filters straddle the redshifted B z band for 
0.5 < z < 0.8. Therefore, we transform the observed 
surface brightnesses to rest-frame B. The procedure is 
described in van Dokkum & Franx (1996). It involves an 
interpolation between F606W and F814W, and is differ- 
ent from applying a "iVcorrection" . For z = 0.583 we 
find 

fiB z = MF814M/ + 0.42(F6061V - F8UW) + 0.84, (6) 
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Fig. 2. — 4" x 4" images (upper) and residuals (lower) after r ' profile fitting for the galaxies in MS 2053-04. Masked 
regions are indicated in black. The postage stamps prove that our results for the early-type galaxies are not suffering from 
misclassifications or bad profile fitting. 



and for z = 0.83 
MSz = [i F814W + 0.01(F606IU - F8UW) + 1.13. (7) 

The F6Q6W - F8UW colors in © and are ob- 
tained using SExtractor (Bertin & Arnouts 1996) with 
fixed apertures of Of! 7 diameter. Extinction corrections 
for both fixed aperture colors and surface brightnesses 
were derived from Schlegel et al. (1998). Another cor- 
rection compensates for cosmological dimming oc (1+z) 4 . 
The final samples are summarized in Table |21 Included 
are the coordinates with respect to the BCG, aperture 
corrected a, r e , fiB, corrected for galactic extinction and 
cosmological dimming, total F814W T magnitudes and 
morphological classifications. 

4. THE FUNDAMENTAL PLANE 



In this section we discuss the FP relation in MS 2053- 
04 and MS 1054-03. For determination of the zeropoint 
and scatter we adopt the slope (a, j3) — (1.20,-0.83) 
that JFK96 found for the local B-band FP. In Fig. 
we show the edge-on view of the FP for both clusters. 
Different symbols indicate different morphological types. 
For comparison the Coma FP is drawn as well. The 
galaxies in MS 2053-04 and MS 1054-03 follow a similar 
FP, but with an offset with respect to Coma. 

4.1. Zeropoint and scatter 

MS 2053-04 

Using the biweight mean (Beers et al. 1990) we 
fit a zeropoint of the FP to the early-type galaxies 
in MS 2053-04. Under the assumption of homology, 
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Fig. 3. — 4" x 4" images (upper) and residuals (lower) after r 1 ' 4 profile fitting for the galaxies in MS 1054-03. Masked 
regions are indicated in black. The postage stamps prove that our results for the early-type galaxies are not suffering from 
misclassifications or bad profile fitting. 
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Visual classification 

Fig. 4. — Visual classification histogram for different best fit- 
ting Sersic numbers. Early-type morphologies generally cor- 
respond to Sersic numbers 3 and 4. A different greyscale is 
used for MS 2053-04 and MS 1054-03. 



the zeropoint shift of the FP traces the mean evolu- 
tion of the galaxy M/L ratio. The observed zeropoint 
offset is Alog(M/L s ) = -0.365 ± 0.037, larger than 
Alog(M/L B ) = -0.280 ± 0.036 found by K97 based on 
older data for a sample of 5 bright galaxies. We will 



return to this issue in Sect. [5] 

We find a biweight scatter for the early-type popu- 
lation in MS 2053-04 as large as 0.134 ± 0.034 in logr e , 
with the error derived from bootstrapping. This is signif- 
icantly larger than the observed scatter of 0.071 around 
the B-band FP of local clusters (JFK96). After subtrac- 
tion of the measurement uncertainties in quadrature, we 
find the intrinsic scatter for the early-type galaxies to 
be 0.124 ± 0.037. We conclude that measurement uncer- 
tainties cannot account for all of the enhanced scatter. 
Not only is the scatter larger than in the local universe, 
it also exceeds the previous result of 0.058 ± 0.018 ob- 
tained by K97. Our new measurements for 4 early-type 
galaxies that were also in the K97 sample give a scat- 
ter of 0.116 ± 0.035 (as opposed to 0.050 ± 0.018 for the 
original K97 data on these 4 galaxies). A larger sam- 
ple and new data on previously studied objects leads to 
the conclusion that the early-type galaxies in MS 2053- 
04 show a considerably larger spread around the FP than 
early-type galaxies in the local universe. 

We next analyzed the scatter of the bulge-dominated 
systems selected by Sersic index (3 and 4), which is a 
more objective method of classifying. The scatter drops 
by about 20% with respect to the classification by eye, 
to 0.111 ± 0.024. Using both visual and quantitative 
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classifications, we find that bulge-dominated systems in 
MS 2053-04 are less tightly confined to a plane than in 
the local universe by a factor 1.5 to 2. 

Zeropoint shifts and scatters from biweight statistics 
are summarized in Tabled Median zeropoint shifts are 
given for comparison. 

MS 1054-03 

For the early-type galaxies in MS 1054-03 the zero- 
point offset agrees well with the previous result from 
vD98 (Alog(M/L B ) = -0.405 ± 0.037 for the new sam- 
ple and A\og(M/L B ) = -0.393 ± 0.040 from vD98). 
The observed scatter is 0.106 ± 0.024, and the intrin- 
sic scatter is 0.086 ± 0.028, consistent with the scatter 
in local clusters (JFK96). Using new measurements of 5 
MS 1054-03 early-type galaxies studied by vD98, we ob- 
tain a biweight scatter of 0.062 ± 0.018, consistent with 
0.047± 0.024 for the original vD98 data on these objects. 
As for MS 2053-04, the scatter decreases by roughly 20% 
if we select bulge-dominated systems by Sersic index (3 
and 4) instead of by eye. 

5. CORRELATIONS WITH OTHER PARAMETERS 

One of the striking effects we found in Sect. I4.ll is the 
large FP scatter for MS 2053-04. Here we investigate the 
cause of this enhanced scatter. In the case of a stellar 
population effect such as variations in age, metallicities 
or dust content, we expect the residual from the FP to 
correlate with color and linestrength indices. We also 
discuss the residual from the FP as a function of envi- 
ronment and investigate the dependence on galaxy mass. 
Hereafter, we again adopt the JFK96 slope and we ex- 
press residuals from the FP as deviations in log(M/ Lb). 
A positive residual means the galaxy has a higher M/L 
than the FP prediction based on its r e and a. The ap- 
proach of measuring offsets along the surface brightness 



axis is physically intuitive since -if we ignore mergers- a 
galaxy only moves along this axis during its lifetime. For 
consistency with our FP analysis, we adopt locally deter- 
mined slopes for all other considered scaling relations as 
well. In Fig.[S]the residual from the FP is plotted against 
various properties for both clusters. Different symbols re- 
fer to different morphological classes. For each panel, the 
probability that a random sample has the same Spear- 
man rank order correlation coefficient as the early-type 
galaxies in our sample, is printed in the corner. 

5.1. The color-magnitude relation 

The most straightforward interpretation for the large 
scatter is a stellar population effect. For age, metallicity, 
and dust trends, we expect lower M/L ratios than the FP 
prediction to correlate with bluer colors than the CMR 
prediction. We refer to Tran et al. (2003) for the color- 
magnitude relation of MS 2053-04 and to van Dokkum 
et al. (2000) for the MS 1054-03 CMR. Similar to our FP 
analysis, we assume that the slope of the relation does 
not evolve and adopt the slope measured in the Coma 
cluster by Bower, Lucey,& Ellis (1992b). After conver- 
sion of the (F606W,F814W) photometry to Johnson U 
and V (using the same procedure as K2000), we fit a 
CMR zeropoint to the early-type galaxies in our sam- 
ples. A positive residual ACMR in Fig. and Fig. 
EJd corresponds to a redder color of the galaxy than the 
CMR prediction based on its total V magnitude. The 
Spearman rank order correlation coefficient points in the 
expected direction for a stellar population effect. How- 
ever, the level of significance is insufficient to confirm 
that galaxies with lower M/L than the FP prediction 
are bluer and the higher (more evolved) ones are redder 
than the CMR prediction. 

5.2. Hj3 linestrength 
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Table 1. MS 2053-04 and MS 1054-03 zeropoint and scatter around FP. 
Earlier results from K97 and vD98 are also in this Table. 



Sample 


# objects 


A log(M/ Lb) 
(biwcight mean) 


Alog(M/L s ) 
(median) 


Scatter in log r e 
(biweight) 


MS 2053-04 early-type 
MS 2053-04 Sersic34 
MS 2053-04 K97 


19 
23 
5 


-0.365 ± 0.037 
-0.382 ± 0.028 
-0.280 ± 0.036 


-0.404 ± 0.037 
-0.404 ± 0.028 
-0.257 ± 0.036 


0.134 ±0.034 
0.111 ±0.024 
0.058 ± 0.018 


MS 1054-03 early-type 
MS 1054-03 Sersic34 
MS 1054-03 vD98 


12 
15 
6 


-0.405 ± 0.037 
-0.368 ± 0.027 
-0.393 ± 0.040 


-0.418 ±0.037 
-0.392 ± 0.027 
-0.405 ± 0.040 


0.106 ±0.023 
0.086 ± 0.026 
0.049 ±0.018 
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Fig. 6. — Residual from the fundamental plane A \og(M/ Lg) plotted against residual from the color-magnitude relation ACMR, residual 
from the H(3 — a relation, distance from BCG and galaxy mass. Different symbols refer to different morphological types. Error bars for 
colors and angular distances are assumed to be smaller than the symbol sizes. For the early-type galaxies, the p-values for statistical 
significance from the Spearman rank order correlation test are printed in the lower right corners. In panels g) and h) the magnitude limit 
J w 21.15 at which serious incompleteness due to uncertainties in the a measurements sets in, is indicated with the solid line. 



8 



Apart from M/L ratio (residual from FP) and color 
(residual from CMR), strengths of absorption lines are 
a valuable tool for tracing stellar population effects. In 
this section we discuss the Hf3 index (expressed in A, 
see Trager et al. 1998). Hf3 is an age-sensitive param- 
eter with only a minor contribution due to metallicity. 
For spectral reduction and derivation of linestrengths we 
refer to Kelson et al. (2003, in preparation). 

The HP — a relation for early-type galaxies in the 
Coma cluster was derived from the H [3q — a relation pre- 
sented by J0rgensen (1999). The HPq index is related 
to the Lick/IDS Hf3 index as Hf3 G = 0.866iJ/3 + 0.485 
(J0rgensen 1997). Assuming a non-evolving slope of the 
HP — a scaling relation, we fit a zeropoint to the rela- 
tion in MS 2053-04 and MS 1054-03. The residuals from 
this relation (positive means stronger HP) are plotted 
against the FP residual in Fig. Et and Fig. Eh- Only 
for MS 2053-04 are the error bars small enough to draw 
robust conclusions. A correlation is present, with con- 
fidence level 95.8%. HP absorption lines are stronger 
for younger stellar populations, and the correlation with 
Alog(M/L_s) could confirm that the scatter in the FP 
is not random noise, but determined by age variations 
among the galaxies. 

5.3. Location in the cluster 

Here we consider if the enhanced scatter reported in 
Sect. 14.11 is related to environment. Clusters of galax- 
ies are not isolated systems, as infall of galaxies from 
the field occurs. FP studies of field early- type galax- 
ies indicate that their stellar populations may be some- 
what younger than those of their counterparts in clus- 
ters (van Dokkum et al. 2001; Treu et al. 2002; van de 
Ven, van Dokkum, & Franx 2003; Rusin et al. 2003; van 
Dokkum & Ellis 2003; Gebhardt et al. 2003). In the 
cluster CL1358+62 van Dokkum et al. (1998) reported 
evidence for disky galaxies to be systematically bluer and 
to show a larger scatter in the CMR at larger radii (with 
the sample ranging to ~ IMpc). As earlier studies of 
both clusters (K97; vD98) were based on one pointing, 
the larger scatter we find could be explained if residu- 
als from the FP increase with distance from the BCG. 
Therefore, we might expect to see a gradient in age and 
larger scatter around the FP going to a larger range in 
cluster radii. For both clusters our samples extend to 
roughly IMpc from the BCG. Fig. 03 and Fig.EF do not 
show a significant correlation. 

5.4. Galaxy mass and selection effects 

Finally we try to explain the range of FP residuals as 
a function of galaxy mass. The mass M in solar units of 
a galaxy is calculated as follows (see JFK96): 



be rewritten as 



log M = 2 logo- + log7- e + 6.07 



(8) 



Fig. Eh and Fig. Eh show residual from the FP against 
logM. Only if we were to include the 6 galaxies 
with a < 100 km s , the trend of lower mass galax- 
ies to have lower M/L ratios than predicted by the 
FP with JFK96 coefficients is significant at the 95% level. 

For a good understanding of Fig. Ek and Fig. Eh, 
we need to take into account that our FP samples are 
magnitude-limited, and not mass-limited. The FP can 



M/L ~ M 



0.28 r -0.07 



(9) 



In the following, we ignore the dependence on r e . Hence 
the residual is given by 



AM/L 



M b s /L bs 

(M/L) 



M° ob f 



FP 



L Gl 



(10) 



For a fixed luminosity we expect all galaxies to fall on 
a line in a plot of A log M/L versus logM. In Fig. 
Ek and Fig. Eh this line is drawn for 1=21.15. At this 
magnitude serious incompleteness due to uncertainties 
in the a measurements sets in. The lowest mass galaxies 
in Fig. Eh and Fig. Eh lie close to the line representing 
the magnitude limit. Low mass galaxies that lie on or 
above the FP would be too faint to allow for accurate 
dispersion measurements and would not enter the FP 
samples. Hence the few remaining low mass galaxies in 
our samples are brighter than their FP prediction based 
on r e and a. 

Selection effects are clearly less relevant at the high 
mass end, and hence we determine the offset and scat- 
ter of the MS 2053-04 FP seperately for the subsample 
of early-type galaxies with M > 10 11 M Q (following van 
Dokkum & Stanford 2003). For this subsample of 8 early- 
type galaxies we derive a zeropoint shift with respect to 
Coma of A\og(M/L B ) = -0.288 ± 0.056 and a scatter 
of 0.132 ± 0.039 in logr e . As expected, this is slightly 
different from the zeropoint shift of A\og(M/ Lb) — 
—0.365 ± 0.037 for the full early-type sample. The ze- 
ropoint offset is now in good agreement with K97, but 
the scatter is still larger. 

If we apply a mass cut for the MS 1054-03 early- 
type galaxies of M > 10 115 Mq, we find a shift in ze- 
ropoint of A \og(M / Lb) = -0.311 ± 0.051 compared 
to Alog(M/is) = -0.405 ± 0.037 for the full early- 
type sample. The scatter of the high mass subsample 
(0.104± 0.030 in logr e ) is similar to that of the full early- 
type sample (0.106 ± 0.023 in logr e ). 

5.5. Summary of correlations 

We did not find evidence that FP residuals are related 
to environment. Instead, stellar population effects are 
playing a role in shaping the FP and selection effects 
in our magnitude-limited samples need to be taken into 
account. 

First, early- type galaxies with stronger HP absorption 
also tend to have lower M/L ratios than predicted by 
their r e and a. This correlation supports the interpreta- 
tion of FP scatter as a measure of age variation among 
early- type galaxies. The larger FP scatter in MS 2053- 
04 therefore reflects a larger spread in relative ages than 
in the local universe. Comparison of FP residuals with 
residuals from the CMR cannot confirm or rule out the 
presence of such a stellar population effect. 

Second, the fact that we do not see galaxies with low 
masses that lie on or above the FP, may be entirely ex- 
plained by selection effects. Only low mass galaxies that 
are brightened by some recent star formation enter the 
FP sample. Their older counterparts are fainter than the 
magnitude limit for the dispersion measurements. 
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Fig. 7. — Evolution of the M/L ratio with redshift. The cross symbols are results from the literature, namely Coma at z = 0.023 (JFK96), 
CL1358+62 at z = 0.33 (K2000), CL0024+16 at z = 0.39 (van Dokkum & Franx 1996), MS 2053-04 at z = 0.583 (K97), MS 1054-03 at 
z = 0.83 (vD98) and J0848+4453 at z = 1.27 (van Dokkum & Stanford 2003). Boxes denote the shift in M/L ratio based on our larger 
samples for MS 2053-04 and MS 1054-03. Single burst models for zj orm = 2 and oo, assuming a Salpeter IMF (x = 2.35) and a range 
of metallicities, are drawn with solid and dotted curves resp. Panel (a) is based on all early-type galaxies, panel (b) shows results for 
the subsample of massive early-type galaxies (M > 10 11 Mq for clusters up to z = 0.583 and M > 10 11,5 Mq for the two higher redshift 
clusters). 



6. EVOLUTION OF M/L RATIO 

Correlations with residual from the Hp — a relation 
show that differences in FP zeropoint can be explained 
by age differences of the stellar population. In this sec- 
tion we study the evolution of the M/L ratio as a func- 
tion of redshift and use this to estimate the mean for- 
mation redshift of cluster early- type galaxies. In Fig. 
[7^, we show the evolution of the M/L ratio with re- 
spect to Coma. Crosses refer to measured zeropoint 
shifts from the literature: Coma at z = 0.023 (JFK96), 
CL1358+62 at z = 0.33 (K2000), CL0024+16 at z = 0.39 
(van Dokkum,& Franx 1996), MS 2053-04 at z = 0.583 
(K97), MS 1054-03 at z = 0.83 (vD98) and J0848+4453 
at z = 1.27 (van Dokkum,& Stanford 2003). The new 
results for MS 2053-04 and MS 1054-03, based on our 
larger samples, are plotted with boxes. Single burst evo- 
lutionary models for a formation redshift Zf orrn = 2 and 
Zform = oo are drawn with a solid and dotted curve re- 
spectively. They represent a galaxy that is fixed in mass 
and whose luminosity evolves as: 

L(t) ~ l/(t-t form ) K (11) 

Here tt orm represents the age of the universe at the 
moment the stars were formed, k depends on the slope 
of the IMF, passband and metallicity. For a Salpeter 
(1955) IMF and -0.5 < [Fe/H] < 0.5 models of Bruzual 
k Chariot (1998), Vazdekis et al. (1996) and Worthey 
(1998) give 0.86 < k b < 1.00. Using the new offsets for 
the two higher z clusters, the single burst model for a 
Salpeter IMF and solar metallicity favoured by the least 
square method has Zf orm ~ 2.261q' 2 o- The la confi- 



dence level was derived from the difference in \ 2 between 
the model and the overall minimum, A\ 2 = X 2 ~ Xmin> 
to which Gaussian confidence levels were assigned (e.g. 
Press et al. 1992). As we showed in Sect. 15.41 the point 
of MS 2053-04 deviates from the earlier result (K97) 
since at low mass only galaxies with low M/L enter the 
magnitude-limited sample. 

To avoid this bias, we apply a mass cut of M > IO^Mq 
to all clusters up to z — 0.583. A mass cut of M > 
10 115 M Q was applied to MS 1054-03 and J0848+4453 
since at these higher redshifts the selection effect sets in 
at a higher galaxy mass. For J0848+4453 two galaxies 
are left after omitting the one low mass outlier. As the 
biweight location estimator is robust against outliers, the 
zeropoint shift only changes slightly for this cluster. We 
obtain the evolution of the M/L ratio as presented in 
Fig. 03. Now the zeropoint of the MS 2053-04 FP fol- 
lows the trend seen for the other clusters in Fig. If 
we constrain the analysis of the evolution in M/L ratio 
to massive early- type galaxies, a simple linear fit gives 
log M/ Lb ~ — 0.47z, agreeing well with earlier determi- 
nations based on smaller samples or lower redshift (see, 
e.g. vD98). The formation redshift favoured by a least 
squares method is Zf orm ~ 2.95^o 46) slightly higher than 
the mean formation redshift for all early-type galaxies. A 
similar formation epoch for early-type galaxies in clusters 
was found by Kelson et al. (2001) based on the evolu- 
tion of Balmer absorption- line strengths with redshift. 
It is remarkable how earlier studies of clusters based on 
smaller samples (see, e.g. K97; vD98) gave similar re- 
sults. 

7. SUMMARY 
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We used visual and quantitative classifications to se- 
lect bulge-dominated systems in our MS 2053-04 and 
MS 1054-03 samples. For MS 2053-04 we find a zero- 
point offset with respect to Coma of Alog(Af/Ls) = 
—0.365 ± 0.037, larger than determined earlier on the 
basis of a smaller sample (K97). The scatter around 
the MS 2053-04 FP is 0.134 ± 0.034 in logr e , enhanced 
with respect to both K97 and the scatter in local clus- 
ters. The FP zeropoint of MS 1054-03 (Alog(M/i B ) = 
—0.405 ± 0.037) agrees well with the earlier result of 
Alog(M/L B ) = -0.393 ± 0.040 by vD98. The scatter 
of 0.106 ± 0.024 in MS 1054-03 is larger than reported 
by vD98 for a smaller sample of MS 1054-03 early-type 
galaxies. However, taking into account measurement un- 
certainties, the scatter is consistent with that in local 
clusters (JFK96). Late-type galaxies also follow the FP 
scaling relation and show a similar scatter around the 
FP as the early-type galaxies. Adding the late-type 
galaxies to the early-type sample results in a scatter 
of 0.136 ± 0.029 for MS 2053-04 and 0.117 ± 0.031 for 
MS 1054-03. The larger samples presented in this paper 
allow us to study correlations with other properties of the 
early- type galaxies. We do not find evidence that the 
formation history depends on environment in the clus- 
ter. No significant correlation of FP residuals with CMR 
residuals was found. The presence of a correlation be- 
tween FP residuals and residuals from the Hf3—a relation 
indicates that stellar population effects are playing a role. 
Assuming non-evolving slopes for all scaling relations, we 
find that galaxies with lower M/L than the FP predic- 
tion tend to show stronger HP than predicted based on 
the HP — a relation. Finally, we show that the lack of 
low mass galaxies on or above the FP may be entirely 
due to selection effects. To avoid a bias induced by the 
magnitude limit of our sample, we focus on the high mass 
end, where selection effects are less relevant. Applying 
a mass cut at M > 10 Mq to all 4 considered clusters 
below z ~ 0.6 and at M > lO n - 5 M to the 2 higher 
rcdshift clusters at z = 0.83 and z — 1.27, increases the 
best fitting formation redshift from Zf orm — 2.26+q |q to 

Zform = 2.95lg;4g. The mass cut at M = lO n Af is well 
below the typical mass of early- type galaxies: galaxies 
with M = 10 11 Mq have dispersions of ~ 168 km s _1 
which is significantly lower than the <7» dispersion of 
early-type galaxies which is 228 ± 14 km s _1 (Kochanek 
1994). 

The implication of this work is that selection effects need 
to be taken into account, especially if the scatter is high. 
The scatter in MS 2053-04 is slightly higher than at low 
redshift; the scatter in the field at high redshift seems to 
be even higher (e.g. Gebhardt et al. 2003; van Dokkum 
& Ellis 2003). Hence those studies are likely to suffer 
from much more significant selection effects. 
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Table 2. FP sample 
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Coordinates with respect to the BCG of MS 2053-04: ID1667 at (20:56:21.4; -04:37:50.8) (J2000). 
b Surfacc brightnesses hb z are corrected for galactic extinction and cosmological dimming. 
Coordinates with respect to the BCG of MS 1054-03: ID4520 at (10:56:59.9; -03:37:37.3) (J2000). 



